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Abstract. We present a possible interpretation for the increase of the angular diameter of the supergiants a Ori (Ml -2 Ia-Ibe) 
and a Her (M5 Ib-II) from the K band to the 1 1 yum region and the high-resolution 1 1 /im spectra without any salient spectral 
features revealed by Weiner et al. (2003a). The angular diameters as well as the high-resolution spectra of a Ori and a Her 
obtained in the 1 1 jim region can be reproduced by a warm water vapor envelope, whose presence in a Ori was revealed by 
Tsuji (2000a) based on the reanalysis of the near-infrared data obtained with the Stratoscope II. While prominent absorption 
due to H 2 can be expected from such a dense, warm water vapor envelope, the absorption lines can be filled in by emission 
from the extended part of the envelope. This effect leads to a significant weakening of the H2O lines in the 1 1 /im region, and 
makes the observed spectra appear to be rather featureless and continuum-like. However, the emission due to H 2 lines from 
the extended envelope leads to an increase of the apparent size in this spectral region. The observed angular diameter and the 
high resolution spectra of a Ori and a Her in the 1 1 /mi region can be best interpreted by the water vapor envelope extending 
to 1.4 - 1 .5 /?* , with a temperature of ~ 2000 K and a column density of H2O of the order of 10 20 cm -2 . 
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1. Introduction 

Recent mid-infrared interferometric observations have revealed 
that the diameters of late-type stars such as M supergiants and 
Mira variables increase from the K band to the N band. Weiner 
et al. (2000, 2003a, hereafter WOO and WHT03a, respectively) 
observed the supergiants a Ori (Ml -2 Ia-Ibe) and a Her (M5 
Ib-II) at 1 1 /mi, using the Infrared Spatial Interferometer (ISI) 
with a narrow spectral bandwidth of 0. 17 cm -1 . They found that 
the uniform disk (UD) diameters of these stars are ~ 30% larger 
than those measured in the K band, which at least for these stars 
are considered to approximately represent the continuum radii, 
as compared to those measured in the optical where scattering 
due to dust particles is more pronounced (see, e.g., discussion 
in Dyck et al. 1992). WOO and WHT03a show that the angular 
diameter of a Ori at 1 1 [im ranges from 53 to 56 mas (43 - 
44 mas in the K band derived by Dyck et al. 1992 and Perrin et 
al. 2004), while that of a Her at 1 1 fim is 39 mas (31-32 mas 
in the K band derived by Benson et al. 1991 and Perrin et al. 
2004). 

Although dust emission is detected for a Ori, the extended 
dust shell is unlikely to be responsible for this increase of the 
diameter for the following reason. The visibility expected for 
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a system consisting of a stellar disk and a very extended dust 
shell is characterized by a steep drop at low spatial frequen- 
cies and a gradual decrease at higher spatial frequencies. The 
amount of the steep drop at low spatial frequencies corresponds 
to the flux contribution of the extended dust shell in the field 
of view, while the stellar disk affects the visibility shape at 
higher spatial frequencies. The inner radius of the dust shell 
around M supergiants such as a Ori is derived to be ~ 50 R+ 
by Danchi et al. (1994). Such a large dust shell is completely 
resolved with the baselines used by WOO and WHT03a (20 - 
56 m), that is, the dust shell does not affect the shape of the 
observed visibility function at these baselines. The presence of 
the dust shell lowers the total visibility at these baselines by 
an amount equal to the fraction of flux coming from the dust 
shell in the field of view. WOO and WHT03a take this effect 
into account in deriving the uniform disk diameters, and show 
that the observed visibilities can be well fitted by a uniform 
disk, but with diameters clearly larger than those measured in 
the K band. For a Her, Danchi et al. (1994) detected a circum- 
stellar dust shell, but more recent spectrophotometric observa- 
tions by Monnier et al. (1998) show that a Her exhibits very 
little dust emission. WHT03a suggest that the dust shell around 
a Her might have evolved since the observations by Danchi et 
al. (1994). Therefore, the increase of the angular diameter ob- 
served in a Ori and a Her cannot simply be attributed to the 
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dust shell, while it cannot be completely ruled out that a pos- 
sible presence of small dust clumps close to the star affects 
the observed visibility, as Bester et al. (1996) suggest based on 
their ISI observations of a Ori. 

WHT03a also present high-resolution spectra of a Ori and 
a Her in the same wavelength range as they selected for the 
ISI observations. These high-resolution spectra were obtained 
with the TEXES instrument mounted on the Infrared Telescope 
Facility with a spectral resolution of ~ 10 5 (Lacy et al. 2002). 
A glance of these TEXES high-resolution spectra (Fig. 2 of 
WHT03a) reveals that neither a Ori nor a Her shows any sig- 
nificant spectral features in the bandpasses used for the ISI ob- 
servations. WHT03a conclude that the interferometric obser- 
vations with ISI were carried out in the bandpasses which ex- 
hibit continuous spectra free from any molecular and/or atomic 
features. That is, the diameters measured in the 1 1 //m region 
which appears to be the continuum are significantly larger than 
the K band diameters, which are also considered to be rather 
close to the diameter measured in the continuum. Therefore, 
the interferometric and spectroscopic observations by WHT03a 
pose a serious problem in interpreting the increase of the an- 
gular diameter from the near-infrared to the mid-infrared. As 
possible explanations for these apparently inconsistent results, 
WHT03a suggest that the density stratification may be very 
nonhydrostatic near the photosphere and/or that the angular di- 
ameters measured in the near-infrared may be affected by a 
presence of hot spots, which may make the apparent size of 
a star smaller than it is. 

Tsuji (2000a, hereafter TOOa) discovered a warm water va- 
por layer around two supergiants a Ori and /i Cep based on 
the reanalysis of the near-infrared spectra obtained with the 
Stratoscope II. The presence of water vapor in a Ori was fur- 
ther confirmed by the detection of weak H2O absorption in the 
6 /mi spectrum obtained with the Infrared Space Observatory 
(ISO) (Tsuji 2000b). TOOa also studied the pure-rotation lines 
of H2O identified in the high-resolution spectrum of a Ori at 
12 fim obtained by Jennings & Sada (1998). Although the ob- 
served equivalent widths of the water lines at 12 /mi are larger 
than those predicted from the photospheric model of TOOa by 
a factor of 2 - 3, the 12 /im water lines are expected to be 
even stronger than observed, if the column density of water 
molecules derived from the near-infrared spectrum is correct. 
TOOa and Tsuji (2000b) point out the possibility that the ab- 
sorption may be filled in by emission from the extended part of 
the water vapor envelope. 

If the warm water vapor envelope is present and the emis- 
sion from it can make the 12 /im H2O lines appear to be weaker 
by the filling-in effect, it may also be the case for the 1 1 /mi 
lines presented in WHT03a. That is, the absorption is filled 
in by the emission resulting from the geometrical extension of 
the warm molecular envelope, making the spectra almost fea- 
tureless. While such featureless mid-infrared spectra may not 
show any hint of the presence of the warm molecular envelope 
and manifest themselves as the continuum, the extended, warm 
molecular envelope can cause the apparent diameter in the mid- 
infrared to be larger than the photospheric diameter. Therefore, 
the warm water vapor envelope may explain the increase of an- 




Fig. 1. The model used in the present work. The photosphere is ap- 
proximated with a blackbody of 7^, with a radius R*. The column 
density in the radial direction (N co {) and the temperature of H 2 gas 
(Tmoi) as well as the radius of the water vapor envelope (R mo \) are the 
input parameters. 



gular diameters from the K band to the N band, simultaneously 
with the spectra observed in the 1 1 /im region. 

In the present paper, we examine this possibility for a Ori 
and a Her. We compare angular diameters as well as spectra 
predicted from a simple model for the warm water vapor enve- 
lope with several observational data: 11 /mi and 12 /im high- 
resolution spectra, the ISO spectrum in the 6 /im region, and 
the angular diameters measured at 1 1 /im. 



2. Modeling of the warm molecular envelope 

Our model consists of the central star, which is represented by 
a blackbody of effective temperature T e s, and a warm molec- 
ular envelope extending to R mo \, as depicted in Fig. 1. As the 
figure shows, the inner radius of the molecular envelope is set 
to be equal to the radius of the star /?*. We adopt an effective 
temperature of 3600 K for a Ori (Tsuji et al. 1994) and 3200 K 
for a Her (Tsuji 1981). The input parameters are the tempera- 
ture of water vapor gas (r mo i), its column density in the radial 
direction (N co {), and the geometrical extension of the envelope 
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(,R mo i in units of /?*). The temperature and the density of the 
water vapor are assumed to be constant in the envelope. 

We first calculate the line opacity due to H2O using the 
HITEMP line list (Rothman 1997), with a Gaussian profile as- 
sumed. We adopt a velocity of 5 km s" 1 for the sum of the 
thermal velocity and the microturbulent velocity. The energy 
level populations of H2O are calculated in local thermodynam- 
ical equilibrium (LTE). The validity of LTE can be examined, 
using order-of-magnitude estimates of collisional and radiative 
de-excitation rates, as adopted by Ryde et al. (2002). The colli- 
sional de-excitation rate C u i is given by C u i ~ Nho~u\v k i, where 
Nn is the density of H atoms, which are assumed to be the 
primary collision partner, <x u i is the cross section, which we 
approximate with the geometrical cross section, and v re i is the 
relative velocity between the H atoms and H2O molecules. As 
we will show below, the column densities of H2O of the warm 
water vapor envelope in a Ori and a Her are derived to be 2 - 
7 x 10 20 cm 2 . The radius of the water vapor envelope is de- 
rived to be 1.4 - 1.5 R+, which is translated into ~ 2 x 10 13 cm 
with a stellar radius of 650 Rq assumed. The number density 
of H2O is then estimated to be 1.0 - 4 x 10 7 cm 4 . The ratio 
of the number density of H atoms to that of H2O molecules ex- 
pected in chemical equilibrium is approximately 10 4 - 10 5 for 
the relevant temperatures and densities. Therefore, the number 
density of H atoms is estimated to be 1 x 10 11 -4 x 10 12 cirT 3 . 
With a geometrical cross section <x u i of 10~ 15 cm 2 and a rel- 
ative velocity v re i of 5 km s -1 assumed, these number densi- 
ties of H atoms lead to collisional de-excitation rates of ~ 50 

- 2000 s -1 . On the other hand, the rate of spontaneous emis- 
sion can be estimated from the Einstein coefficients A u \. For 
the H2O molecule, the ranges of A u \ are approximately ~ 10" 4 

- 3 x 10 2 s _1 , ~ 10~ 5 - 3 x 10 2 s _1 , and ~ 10 s _1 for the 11 

- 12 pm region, K band, and L' band, respectively. Therefore, 
for the wavelength regions that we will discuss below, the as- 
sumption of LTE is valid for weak and moderately strong H2O 
lines, while non-LTE effects may not be negligible for strong 
lines. However, a quantitative assessment of non-LTE effects is 
beyond the scope of the present paper, and we assume LTE for 
the H2O lines considered here. 

Once the line opacity is calculated, the intensity distribution 
at the wavenumber a> can be calculated as 



3. Comparison with the observed spectra and 
angular diameters: a Ori 

Using the above model for the warm H2O envelope, we com- 
pare the synthetic spectra and the visibilities with the observed 
data, which include the 1 1 pm spectra presented in WHT03a, 
the 12 /im spectra obtained by Jennings & Sada (1998), the 
6 pm spectrum in Tsuji (2000b), and the interferometric obser- 
vations at 1 1 pm obtained by WOO and WHT03a. We change 
the input parameters (T mo \, R mo i, 7V co i), and search for the com- 
bination which can simultaneously reproduce these observa- 
tional results. The range of the input parameters used in the 
calculations is as follows: r mo i (K) = 2150, 2100, 2050, 2000, 
1950, 1900, 1800, 1700, 1600, 1400, and 1200, N col (cnr 2 ) = 
1 x 10 19 , 1 x 10 20 , 2 x 10 20 , 5 x 10 20 , 1 x 10 21 , and 1 x 10 22 , 
flmoi (R-k) = 1.3, 1.4, 1.45, 1.5, 1.55, 1.6, and 1.7. It has turned 
out that the aforementioned observational data for a Ori can 
be best reproduced with r mo i = 2050 K, R mo \ = 1.45 /?*, 
A^ co i = 2 x 10 20 crrT 2 . The uncertainties of these values are 
estimated by changing the parameters around the best-fit pa- 
rameter set by small amounts. We estimate that the uncertain- 
ties of the gas temperature, the radius, and the column density 
of H2O are ±100 K, ±0.1 /?*, and a factor of ~ 2, respectively. 
The gas temperature and the column density of H2O molecules 
derived in the present work are in rough agreement with those 
derived by TOOa. In the following subsections, we will discuss 
the comparison for each observational data set. 

3.1. 11 /urn spectra 

We calculate synthetic spectra in the 11 pm region with a 
wavenumber interval of 0.001 cm"', and then convolve with a 
Gaussian profile which represents the effects of the instrument 
and the macroturbulent velocity. The spectral resolution of the 
TEXES instrument is 10 5 , which translates into an instrumen- 
tal broadening of 3 km s -1 . Jennings & Sada (1998) derived a 
macroturbulent velocity of 12 km s -1 for a Ori. Thus, the syn- 
thetic spectra are convolved with a Gaussian with a FWHM of 
V3 2 + 12 2 = 12.4 km s _1 , which corresponds to 0.037 cirT 1 at 
897 cirT 1 . In order to take the continuous dust emission from 
the circumstellar dust shell into account, the convolved and nor- 
malized spectrum F w is diluted as follows: 



I(p, 10) = BoX^moiXl - e~ T ") + BJT eff )e- T " circ(p/R*), (1) fluted = (1 _ /dust)jFtj + /di 



(3) 



where p is the impact parameter, t m is the optical depth along 
the ray at p, and circ(p/^*) takes a value of 1 for p < R* and 
elsewhere (see Fig. 1). In order to obtain the spectrum, the 
emergent flux at each wavenumber, F u , is calculated by 



= 2n \ 

Jo 



I(p, <x>) p dp, 



(2) 



Vl - (P/Rmoi) 2 With G 



where p is defined as p — cos 6 
defined in Fig. 1. With an appropriate wavenumber interval, 
we calculate an intensity profile and a flux at each wavenum- 
ber. The corresponding (monochromatic) visibility is calcu- 
lated from this intensity profile using the Hankel transform 
(two-dimensional Fourier transform for axisymmetric objects). 



where /fluted is the final spectrum, and /dust is the fraction of 
the flux contribution of the circumstellar dust shell. WHT03a 
derived /d USt ~ 0.44 for a Ori at 1 1 pm, and we adopt this value. 
Figure 2 shows a comparison between the synthetic 1 1 pm 
spectra from the best-fit model for a Ori and the observed spec- 
tra presented in WHT03a. The dots in the figure represent the 
observed data, which were read off Fig. 2 in WHT03a, while 
the solid lines represent the synthetic spectra, which are red- 
shifted by 39 km s -1 to match the observations (see WHT03a). 
Since our model is too simple to deal with the complicated 
physical processes in the outer atmosphere, we cannot expect 
very good quantitative agreement for individual spectral fea- 
tures. However, Fig. 2 illustrates that the model calculation can 
reproduce the absence of conspicuous spectral features as well 
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Fig. 2. Spectra in the 1 1 pm region. The dots represent the spectra of 
a Ori presented in WHT03a, while the solid lines represent the cal- 
culated spectra from the best-fit model with T mo \ = 2050 K, R mo \ = 
1.45 R+, and A^oi - 2x 10 20 cirr 2 . The synthetic spectra are convolved 
with a Gaussian with a FWHM of 0.037 cirr 1 to account for the ef- 
fects of the instrument as well as of the macroturbulent velocity, and 
redshifted by 39 km s -1 to match the observations. The positions of 
the H 2 lines whose intensity at 2000 K is stronger than 3 x 10~ 23 cm 
molecule -1 are marked with ticks. These line positions are taken from 
the HITEMP database, and are also redshifted by 39 km s with re- 
spect to the rest wavenumber. 



as the depths of the small features. Although a number of pure- 
rotation lines of H2O are present as shown with the ticks in the 
figure, they appear neither as absorption nor as emission. The 
reason for the absence of salient spectral features is that the ab- 
sorption lines are filled in by the emission from the outer part of 
the envelope. In order to illustrate this effect, we show synthetic 
spectra at the disk center (p — 1) and near the limb of the warm 
H2O envelope (p = 0.3) in Fig. 3 (top panel), together with 
the resulting spectrum (emergent flux) expected for the whole 
system (middle panel), and the optical depth due to H2O in the 
radial direction (bottom panel). As the bottom panel shows, the 
optical depth due to H2O ranges mostly from ~ 0.1 to ~ 1 with 
several peaks reaching ~ 10, which means that the warm water 
vapor envelope is not optically thick. In such a case, the ob- 
server see the water vapor envelope with the hotter star in the 
background in a line of sight with p < R+, which leads to ap- 








Fig. 3. Top panel: Spectra predicted from the best-fit model for a Ori 
at the stellar disk center (p. = 1, upper spectrum) and near the limb 
of the warm H 2 envelope (p = 0.3, lower spectrum) in the region 
around 11.1494 pm. The continuum level is shown with the dashed 
line. The flux contribution of dust emission is not included in the 
spectra shown. Middle panel: The resulting spectrum (emergent flux) 
expected for the whole object (stellar disk + warm water vapor enve- 
lope), calculated from the intensity spectra such as shown in the upper 
panel, using Eq. (2). Bottom panel: Optical depth in the radial direc- 
tion due to H2O lines in the same spectral region. The spectra and the 
flux shown in the top and middle panels as well as the optical depth 
shown in the bottom panel are redshifted by 39 km s" 1 with respect to 
the rest wavenumber to match the observed spectra shown in Fig. 2, 
but not convolved with the Gaussian representing the instrumental ef- 
fect and the macroturbulent velocity. 



pearance of absorption lines. The spectrum at the disk center 
shown in the top panel illustrates this case, where absorption 
features are observed at the positions of the H2O lines. On the 
other hand, in a line of sight with R* < p < R mo i, emission 
lines are observed, as the spectrum near the limb (p. = 0.3) il- 
lustrates. Since the spectrum (i.e. emergent flux) observed for 
the object (central star + water vapor envelope) is obtained by 
integrating the intensity over the area that the object projects 
onto the plane of the sky, the absorption expected from the re- 
gion with p < R* and emission expected from the outer region 
with R* < p < R mo i almost cancel out in the resulting spec- 
trum, leading to the nearly featureless spectra as shown in the 
middle panel. 



3.2. 12 pirn spectra 

Figure 4 shows the synthetic and observed spectra in the 12 pm 
region. The observed spectra, which were read off Fig. 1 in 
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Fig. 4. Spectra in the 12 /mi region. The dots represent the spectra 
of a Ori observed by Jennings & Sada (1998), while the solid lines 
represent the calculated spectra of the best-fit model. The synthetic 
spectra are convolved with a Gaussian with a FWHM of 0. 1 cm" 1 to 
match the resolution of the spectrometer used by Jennings & Sada 
(1998). The ticks represent the positions of H 2 lines whose intensity 
at 2000 K is stronger than 3 x 10~ 23 cm molecule 1 . The observed and 
predicted spectra as well as the H 2 line positions are plotted in the 
laboratory frame. The H 2 lines identified by Jennings & Sada (1998) 
are shown with the arrows. 



Jennings & Sada (1998), are represented with the dots. The 
pure-rotation lines of H2O identified by Jennings & Sada 
(1998) are marked with the arrows in the figure. The synthetic 
spectra, which are represented with the solid lines, were calcu- 
lated with a wavenumber interval of 0.001 cm -1 , and convolved 
with a Gaussian with a FWHM of 0. 1 cm 1 , which corresponds 
to the resolution of the spectrometer used by Jennings & Sada 
(1998). As in the case of the 11 /mi spectra discussed above, 
it is necessary to include the dust emission from the circum- 
stellar dust shell in the synthetic spectra. The flux contribution 
of the dust shell is estimated to be 35% at 12 /mi by Jennings 
& Sada (1998). Therefore, the calculated spectra are diluted 
with /d ust = 0.35. Figure 4 shows that the synthetic spectra can 
reproduce the absence of strong spectral features, and the ob- 
served depths of the weak features are also reproduced to some 
extent. The filling-in effect due to the outer part of the warm 
H2O envelope weakens the absorption lines significantly, as in 




6.5 



Fig. 5. Spectra in the 6 /mi region. The thin solid line represents 
the ISO SWS spectrum of a Ori, while the thick solid line repre- 
sents the synthetic spectrum predicted from the best-fit model. The 
synthetic spectrum is convolved with a Gaussian with a FWHM of 
1.0 cm - '. The ticks mark the H 2 absorption features identified by 
Tsuji (2000b). 



the case of the 1 1 /mi spectra. The absorption features observed 
at 815.4 cirT 1 and 815.95 cirT 1 are affected by the blend of 
OH lines, which are not included in the present calculation of 
the synthetic spectra. It should also be noted that TOOa shows 
that the 12 /mi H2O lines may partially originate in the photo- 
sphere. If the detailed photospheric structure were incorporated 
in our model, the H2O absorption lines would be stronger than 
shown in Fig. 4, which might improve the agreement with the 
observed spectra. 

3.3. 6 /urn ISO spectrum 

Figure 5 shows a comparison between synthetic spectra 
and the spectrum obtained with the ISO Short Wavelength 
Spectrometer (SWS) in the 6-7 /mi region where weak ab- 
sorption due to the H2O V2 fundamental bands was identified 
by Tsuji (2000b). The thin solid line represents the spectrum 
of a Ori (observed on 1997 October 8 UT) retrieved from the 
ISO data archive, while the synthetic spectrum from the best- 
fit model, which is convolved with a Gaussian with a FWHM 
of 1.0 cm -1 to match the resolution of the ISO SWS spec- 
trum (R =2 1600), is shown with the thick solid line. No di- 
lution due to the circumstellar dust emission is assumed in 
this wavelength region. Figure 5 demonstrates that the model 
can fairly reproduce the H2O absorption features observed 
with ISO SWS, which are marked with the ticks in the figure. 
Therefore, we conclude that our warm H2O envelope model is 
also consistent with the 6 /<m H2O spectrum of a Ori observed 
with ISO SWS. 



3.4. Angular diameter in the 1 1 /im region 

The comparison between the synthetic spectra and those ob- 
served in the 11 /mi, 12 /mi, and 6 /mi regions demonstrates 
that the warm water vapor envelope extending to ~ 1.45 /?* 
with a temperature of 2050 K and an H2O column density of 
2 x 10 20 cirT 2 can reproduce the observed spectra. The calcu- 
lation of the synthetic spectra at 1 1 /mi shows that the pres- 
ence of the warm water vapor envelope can yield a featureless, 
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Fig. 6. a: Visibility calculated from the best-fit model for a Ori in 
the region around 11.1494 /im. The visibility shown here is calcu- 
lated for a projected baseline length of 30 m. b: Uniform disk di- 
ameter calculated from the spectrally convolved visibility shown in a 
for a projected baseline length of 30 m. In a and b, two bandpasses 
used by WHT03a are marked with the dashed lines. The bandpasses 
used for the other observations of a Ori lie between those shown (see 
WHT03a). The range of the diameters measured in these bandpasses 
is shown as the hatched region. The dotted line represents the pho- 
tospheric angular diameter adopted in the calculation. A flux contri- 
bution of 44% from the dust shell is included in the calculations of 
the model visibility, but not in the computation of the uniform disk 
diameter. Both plots (a and b) are redshifted by 0.117 cirr' with re- 
spect to the rest wavenumber, which corresponds to a radial velocity 
of 39 km s , to match the observed spectra shown in Fig. 2. c: The 
normalized intensity profile at 896.85 cm" 1 , which is spectrally con- 
volved with the same top-hat function as used in the calculation of 
the spectrally convolved visibility shown in a. d: The spectrally con- 
volved visibility squared at 896.85 cirr 1 (solid line), plotted together 
with the visibility points observed by WHT03a (filled diamonds), e: 
The same spectrally convolved visibility squared as in d, but at higher 
spatial frequencies. 



continuum-like spectrum in the bandpasses used in the ISI ob- 
servations by WHT03a. If the featureless spectra observed at 
1 1 //m are a result of the filling-in effect due to the emission 
from the outer part of the warm water vapor envelope, the an- 



gular diameter in this wavelength region can appear larger than 
the photospheric diameter. This possibility can be examined 
by computing the visibility from the model intensity profile at 
each wavelength. 

We calculate the intensity profile at a wavenumber inter- 
val of 0.001 cm" , and then from this monochromatic inten- 
sity profile, the monochromatic visibility is obtained using the 
Hankel transform. The monochromatic visibility is convolved 
with an appropriate response function which represents the 
spectral resolution of the ISI observations by WHT03a. The 
spectrally convolved visibility is calculated as follows: 



VconvC/, 10) 



— I V(f, aj')S(a)' - a>) dcJ , 



(4) 



where V(f, a>) is the monochromatic visibility at the spatial fre- 
quency / and the wavenumber a>, V cow/ (f, a>) is the spectrally 
convolved visibility, and S(a>) is the spectral response func- 
tion of the bandpass. In the calculation discussed here, S(a>) is 
assumed to be a top-hat function with a width of 0.17 cm , 
which is the bandwidth used for the ISI observations. We as- 
sume that the Zf-band angular diameter represents the photo- 
spheric diameter. Dyck et al. (1992) and Perrin et al. (2004) 
derived the Zf-band uniform disk diameter of 43 - 44 mas for 
a Ori, and we adopt 44 mas as the photospheric diameter. Once 
the visibility function resulting from the stellar disk and the 
warm H2O envelope (not including the dust shell) is calculated, 
the uniform disk diameter at each wavenumber can be derived 
for a given baseline length. We compute the uniform disk di- 
ameter for a baseline length of 30 m, which is the mean of the 
baseline lengths used in the ISI observations by WHT03a. 

Figures 6a and 6b show the calculated visibility and the 
uniform disk diameter of the best-fit model for a Ori in the re- 
gion around 1 1 . 1494 /im. As mentioned in Sect. 1 , the presence 
of the extended dust shell lowers the visibility by an amount 
equal to the fraction of the flux contribution of the dust shell. 
Therefore, the visibility resulting from the stellar disk and the 
warm H2O envelope is lowered by a factor of 0.56, which ac- 
counts for the flux contribution of the stellar disk in this wave- 
length region. Note, however, that the uniform disk diameter 
is computed from the visibility excluding the dust shell, be- 
cause the effect of the presence of the dust shell is already taken 
into account in the determination of the uniform disk diameters 
by WOO and WHT03a. Therefore, the uniform disk diameter 
shown in Fig. 6 can readily be compared with those observa- 
tionally derived by WOO and WHT03a. As Fig. 6b illustrates, 
the uniform disk diameter in this spectral region is larger than 
the photospheric diameter of 44 mas (dotted line in the fig- 
ure), and the predicted diameter in the ISI bandpasses between 
896.7 and 897.0 cirr 1 is in agreement with the result obtained 
by WHT03a. 

Figure 6c shows the intensity profile at the center of the 
range of the bandpasses used by WHT03a, while Fig. 6d shows 
the corresponding visibility squared as a function of spatial fre- 
quency, together with the observed values presented in Fig. 1 
of WHT03a. Fig. 6d demonstrates that the predicted visibility 
squared is somewhat lower that those observed. The observed 
visibility squared can be well fitted with a uniform disk diam- 
eter of 52.66 mas (WHT03a), while the predicted uniform disk 
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diameter is 56.6 mas. However, Weiner et al. (2003b) show that 
the uniform disk diameter measured for a Ori fluctuates with 
an amplitude of ± ~ 1.5 mas, and given this temporal variation 
of the angular diameter of a Ori on the one hand and the sim- 
plicity of our model on the other hand, the agreement between 
the observed and predicted visibility squared can be regarded 
as fair. Figure 6e shows the same predicted visibility squared 
at higher spatial frequencies, which correspond to a projected 
baseline as long as ~ 100 m. Observations with such long base- 
lines may be realized by ISI in the near future, and the visibility 
shape at such high spatial frequencies will be useful for further 
examining the model for the warm water vapor envelope. 

We perform the same calculation for the other two spec- 
tral regions observed by WHT03a. Figure 7 shows the pre- 
dicted uniform disk diameter in the regions around 1 1 .0856 /mi 
and 11.1713 /im. In the 1 1.0856 /im region, the predicted uni- 
form disk diameter ranges from 48 to 53 mas within the band- 
pass used by WHT03a, which are marked with the dashed 
lines in the figure. These predicted diameters are systemati- 
cally lower than the observed value of 54.14 ± 0.52 mas. Given 
the above mentioned temporal variation of the angular diame- 
ter of a Ori, however, this slight discrepancy is not regarded as 
serious disagreement. The predicted diameter in the ISI band- 
pass in the 11.1713 //m region ranges from 51 to 56 mas (lower 
panel), which is in good agreement with the observed value of 
54.20 ± 0.46 mas. Thus, our simple model of the warm H2O 
envelope can reasonably explain the observed increase of the 
angular diameter of a Ori and the spectra in the near-infrared 
(6-7 /im) and in the mid-infrared (11 /mi and 12 /im). 

3.5. Uniform disk diameter in the K and V bands 

While the emission from the warm H2O envelope in a Ori can 
explain the mid- and near-infrared spectra as well as the angu- 
lar diameter measured with the narrow bandpasses at 1 1 /mi, a 
question may arise: if such a dense H2O envelope is present, 
can it affect the angular diameter in the near-infrared? TOOa 
demonstrates that the absorption feature at 1.9 /im observed in 
a Ori can be explained by the warm H2O envelope, confirming 
that the presence of the warm molecular envelope can be de- 
tected by detailed stellar spectroscopy in the near-infrared. If 
the angular diameter measured in the near-infrared were also 
significantly affected by the warm molecular envelope, our as- 
sumption that the Zf-band angular diameter of a Ori represents 
the photospheric diameter could not be justified, and hence it 
might not be concluded that the warm H2O envelope is re- 
sponsible for the increase of the apparent diameter from the 
near-infrared to the mid-infrared. Moreover, the L'-band angu- 
lar diameter of a Ori measured by Chagnon et al. (2002) is 42 
- 43 mas, which is very close to the /f-band diameter. This 
suggests that the effect of the warm molecular envelope on the 
L'-band diameter should not be prominent. 

We examine the effect of the warm H2O envelope on the an- 
gular diameter measured in the near-infrared by performing the 
same calculations as in the previous subsection, for the K and 
L' bands. Only H2O lines were included in the calculation, and 
we approximate the K- and L'-band filter response functions 




Fig. 7. Uniform disk diameter around 11.0856 //m (upper panel) and 
1 1.1713 /im (lower panel) calculated from the best-fit model for a Ori. 
Both plots are redshifted by 0. 1 17 cm"' in wavenumber, which corre- 
sponds to a radial velocity of 39 km s _1 , to match the observed spectra 
shown in Fig. 2. The bandpasses used by WHT03a are marked with 
the dashed lines, and the ranges of the angular diameters measured by 
WHT03a are represented as the hatched regions. The dotted lines rep- 
resent the photospheric angular diameter adopted in the calculation. 



with top-hat functions centered at 2.15 //m with AA = 0.5 //m 
and at 3.8 /im with AA = 0.6 //m, respectively. The uniform disk 
diameters predicted in the K and L' bands have turned out to be 
44.5 mas and 45.6 mas, respectively, and these values are very 
close to the 44 mas which we adopted as the photospheric an- 
gular diameter. The predicted angular diameter in the L' band 
is also in agreement with the values derived by Chagnon et al. 
(2002). 

The optical depth due to the H2O lines in the K and L' 
bands ranges from ~ 0.01 to ~ 1 for the gas temperature and 
the H2O column density of the best-fit model (L mo i = 2050 K 
and 7Y co i = 2 x 10 20 cm -2 ), while Fig. 3 shows that the opti- 
cal depth in the 11.1494 region mostly ranges from ~ 0.01 to 
— 1. It means that the H2O envelope is not yet totally optically 
thick, and emission from the star can also be observed. In such 
a case, the effect of emission from the warm H2O envelope on 
the angular diameter is not solely governed by emission from 
the envelope, but by the intensity ratio between the star and the 
envelope. With the emission from the star and the warm H2O 
envelope represented with blackbodies of T e g = 3600 K and 
L mo i = 2050 K, respectively, the intensity ratio is expressed 
as -B / i(L mo i)/B / i(r e ff). This intensity ratio is small in the near- 
infrared, but increases toward longer wavelengths. Therefore, 
the effect of the warm H2O envelope on the angular size is 
much less prominent in the K and L' bands than in the 1 1 //m 
region. In the TT-band spectrum of a Ori, the absorption due to 
CO and CN originating in the photosphere is present. However, 
the effects of these molecular absorption features on the K- 
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Fig. 8. Spectra in the region around 11.1494 jim. The dots repre- 
sent the observed spectrum of a Her, while the solid line represents 
the synthetic spectrum from the best-fit model. The synthetic spec- 
trum is convolved with a Gaussian with a FWHM of 0.037 cm" , and 
blueshifted by 0.084 cirr 1 in wavenumber, which corresponds to a ra- 
dial velocity of -28 km s" 1 (see WHT03a), to match the observation. 
The positions of the H 2 lines whose intensity at 2000 K is stronger 
than 3 x 10~ 23 cm molecule -1 are marked with ticks. These line posi- 
tions are taken from the HITEMP database, and are also blueshifted 
by 28 km s" 1 with respect to the rest wavenumber. 



band angular diameter are expected to be minor, because the 
geometrical thickness of the photosphere of a Ori predicted by 
the classical hydrostatic model is only 9% of the stellar contin- 
uum radius (see TOOa), and is much smaller than the radius of 
the warm H2O envelope of ~ 1 .45 R+. 

4. Comparison with the observed spectra and 
angular diameters: a Her 

For a Her, we compare the 1 1 //m spectrum presented in 
WHT03a and the angular diameter measured in this spectral 
range. With an effective temperature of 3200 K adopted for 
a Her, we search for the combination of the input parameters 
which can best reproduce the spectra and the angular diame- 
ter observed in the 1 1 yum region. The parameters of the best- 
fit model are found to be r mo i = 2000 K, R mo \ = 1.4 R+, and 
Ncoi - 1 x 10 20 crrT 2 . We estimate that the uncertainties of the 
gas temperature, the radius, and the column density of H2O are 
±100 K, ±0.1 R*, and a factor of ~ 2, respectively. 

Figure 8 shows a comparison between the observed spec- 
trum (dots) and the calculated one from the best-fit model (solid 
line). As in the case of a Ori, the synthetic spectrum is con- 
volved with a Gaussian which represents the effects of the in- 
strument as well as of the macroturbulent velocity in the at- 
mosphere of a Her. Tsuji (1988) analyzed the high-resolution 
spectra of the CO first overtone bands of M giants with the use 
of photospheric models, and detected an absorption excess in 
the low excitation CO lines, which he concluded originates in 
the warm molecular envelope. Based on the analysis of this ab- 
sorption excess, Tsuji (1988) found that the turbulent velocity 
in the warm molecular envelope of the M giants in his sam- 
ple, including a Her, can be as large as 10 km s -1 . Such a large 
turbulent velocity is in fact observed in the 12 jjm spectrum of 
a Ori, as mentioned in Sect. 3.1. In the present work, we ten- 
tatively assume the same macroturbulent velocity in a Her as 
in a Ori, that is, 12 km s -1 . Since a Her shows no dust emis- 
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Fig. 9. Visibility and uniform disk diameter calculated from the best-fit 
model for a Her in the region around 1 1.1494 yum. The visibility and 
the uniform disk diameter shown are calculated for a projected base- 
line length of 30 m. No flux contribution of the dust shell is assumed 
in the calculation of the model visibility and uniform disk diameter. 
The dotted line represents the adopted photospheric angular diame- 
ter, while the range of the angular diameter measured by WHT03a is 
shown as the hatched region. Both plots are blueshifted by 0.084 cm" 1 
with respect to the rest wavenumber to match the observation. 



sion feature in the N band (Monnier et al. 1998), no dilution 
effect due to dust emission is included in the calculation for 
a Her. A glance of Fig. 8 reveals that the model can reproduce 
the observed continuum-like spectrum nearly free from salient 
features and, to some extent, the depths of the fine spectral fea- 
tures. 

We calculate the visibility and the uniform disk diameter 
in this spectral region. In the calculation of the visibility, we 
assume that the photospheric angular diameter can be repre- 
sented by the /T-band angular diameter. Benson et al. (1991) 
and Perrin et al. (2004) derived ^T-band angular diameters of 
31-32 mas for a Her, and we adopt an angular diameter of 
31 mas as the photospheric angular diameter in our calcula- 
tion. The uniform disk diameter is derived from the calculated 
visibility for a baseline length of 30 m. Figure 9 shows that 
the uniform disk diameter is much larger than the photospheric 
diameter of 31 mas in the spectral region at issue. In the band- 
pass used in the diameter measurements by WHT03a, which 
are marked with the dashed lines in the figure, the angular di- 
ameter ranges from 38 to 42 mas, which is in agreement with 
the 39.32 ± 1.04 mas measured by WHT03a. 

We also examine the effect of the warm H2O envelope on 
the diameter in the K and L' bands. The uniform disk diam- 
eters in the K and L' bands predicted from the best-fit model 
for a Her are 33.3 mas and 33.9 mas, respectively, and the pre- 
dicted L'-band diameter is in agreement with the observed di- 
ameter of 31.04 ± 0.26 mas derived by Chagnon et al. (2002). 
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This confirms that the TT-band and L'-band diameters are not 
significantly affected by the presence of the warm H2O enve- 
lope, and that the use of the /f-band diameter as the photo- 
spheric diameter is reasonable. It should also be noted that the 
parameters of the best-fit model for a Her are in rough agree- 
ment with those derived by Tsuji (1988) based on the analysis 
of the absorption excess in the low excitation lines of the CO 
first overtone bands. 

5. Discussion 

We have shown that the increase of the angular diameters from 
the near-infrared to 1 1 //m observed in the supergiants a Ori 
and a Her can be explained by the presence of a warm H2O en- 
velope, and that dense H2O gas with a temperature of ~ 2000 K 
extends to 1.4 - 1.5 R+. These results are in good agreement 
with those recently obtained by Perrin et al. (2004), who could 
reproduce the angular diameter of a Ori measured in the K and 
L' bands as well as that at 11.15 /mi, using a spherical gaseous 
envelope model without a line-by-line calculation of molecu- 
lar opacities. They derived the temperature and the radius of 
such a gaseous envelope of a Ori to be 2055 K and 1 .33 /?*, re- 
spectively, with optical depths of 0.06, 0.026, and 2.33 in the K 
band, L' band, and at 11.15 //m, respectively. Since our model 
predicts the optical depth as a function of wavenumber, we av- 
erage the optical depth predicted from our best-fit model for 
a Ori in the K and L' bands. The averaged optical depths of 
the warm water vapor envelope predicted from our model are 
0.05 and 0.026 for the K and L' bands, respectively, which is 
in good agreement with the above values derived by Perrin et 
al. (2004). For the optical depth in the 11.15 //m region, Fig. 3 
shows that the optical depth predicted from our best-fit model 
for a Ori ranges from ~ 0.01 up to 10, which is also in rough 
agreement with the value derived by Perrin et al. (2004). 

From UV observations, M supergiants are known to have 
chromospheres as hot as T e ~ 8000 K. Gilliland & Duprees 
(1996) obtained the first image of a Ori at 2550 A with the 
Hubble Space Telescope, and found that the chromospheric ex- 
tension of a Ori is about 3 times as large as the size measured 
in the K band. However, the VLA observations of a Ori by Lim 
et al. (1998) revealed the presence of cooler gas with tempera- 
tures of 1300 - 3400 K at 2 - 7 R*. Lim et al. (1998) conclude 
that the hot chromosphere and the cool gas coexist, but with the 
latter being the more dominant component. Harper et al. (2001) 
constructed a semiempirical model to explain these VLA ob- 
servations, and their one-dimensional model has a tempera- 
ture distribution which first decreases outward from the pho- 
tosphere, and rises to the maximum of ~ 3800 K at ~ 1 .45 R*, 
and then decreases again. While the temperature and the radius 
of the warm H2O envelope we derived for a Ori (~ 2050 K 
at ~ 1.45 /?*) seem to be in conflict with this semiempirical 
model, it should be noted that Harper et al. (2001) suggest that 
their one-dimensional model might represent the average of an 
inhomogeneous structure, where the cool gas (~ 2000 K) and 
the hot plasma (~ 8000- 10 000 K) coexist at ~ 1.45/?*. Then 
it is implied that the warm H2O envelope which we modeled in 
the present work may be part of the cool component detected 
by the VLA observations and modeled by Harper et al. (2001). 



The increase of the angular diameter is detected not only in 
supergiants but also in Miras and non-Mira M giants. Weiner 
et al. (2003b) as well as WOO and WHT03a observed the Mira 
variables o Cet, R Leo, and^f Cyg with the same observational 
technique as applied to a Ori and a Her, and found that the 
angular diameters of these Mira variables are roughly twice 
as large as those measured in the K band. As in the cases of 
a Ori and a Her, they used bandpasses which appear not to 
be contaminated by H2O or other spectral lines. However, as 
we have shown above, the interpretation of the mid-infrared 
spectra of late-type (super)giants is complicated by flux con- 
tribution from the extended outer atmosphere. The increase of 
diameter is detected not only in the mid-infrared, but also in the 
L' band. Menneson et al. (2002) show that the L'-band angular 
diameters of Mira variables as well as semiregular M giants 
are by a factor of 1.2 - 2.0 larger than the Zf'-band diame- 
ters. The increase of the angular diameters from the TT'-band 
to the L'-band and the 1 1 fim region observed in the Mira vari- 
ables and semiregular M giants may also be explained by the 
warm molecular envelope, whose presence in these classes of 
objects is detected by the analyses of infrared molecular spectra 
(e.g., Tsuji et al. 1997, Yamamura et al. 1999, Cami et al. 2000, 
Matsuura et al. 2002). This possibility will be further studied in 
a forthcoming paper (Ohnaka 2004). 

Although our ad hoc model for the warm molecular enve- 
lope can reproduce the observed spectra as well as the angular 
diameter increase from the near-infrared to the mid-infrared, it 
is not a unique solution. The hypothesis of the warm molecular 
envelope should be further examined by comparing with spec- 
troscopic and interferometric data in other wavelength regions. 
For example, interferometric observations in many more band- 
passes in the 11 //m region would provide a more complete 
picture of the wavelength dependence of the angular diame- 
ter, which would be a further constraint for modeling the warm 
molecular envelope. Interferometry with even higher spectral 
resolution would also be very useful for testing the hypothe- 
sis proposed in the present work. Even if the basic picture of 
the warm water vapor envelope is correct, our simple model is 
not sufficient to understand detailed physical properties of the 
warm molecular envelope such as temperature and density dis- 
tributions, and therefore, the physical parameters we derived 
above should be regarded as representative values of the real 
molecule forming region. 

Furthermore, the physical mechanism responsible for the 
formation of the warm molecular envelope remains to be an- 
swered. In the coolest and the most luminous objects such as 
Mira variables, the levitation of the atmosphere due to stel- 
lar pulsation may lead to a density enhancement in the outer 
atmosphere, and hence creating an environment favorable for 
molecules to form, as Helling & Winters (2001) discuss for 
carbon-rich objects. However, it is not yet clear whether such 
a mechanism can operate in M (super)giants with higher ef- 
fective temperatures and much smaller variability amplitudes. 
For example, Tsuji et al. (1997), Matsuura et al. (1999), and 
Tsuji (2001) detected water vapor in the near-infrared spectra 
of early M giants and a late K giant, whose photospheres had 
been deemed to be too hot for water vapor to form. Ryde et al. 
(2002) detected H2O pure-rotation lines at 1 1 - 12 //m even in 
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the K1.5III giant a Boo. They argue, however, that these H2O 
lines do not originate in the warm molecular envelope, but in 
the outer layers of the photosphere where temperature may de- 
viate from that predicted by classical, hydrostatic photospheric 
models, although the mechanism responsible for such a devia- 
tion remains ambiguous. In any case, the physical properties of 
the region between the upper photosphere and the expanding, 
cold circumstellar envelope have been increasingly probed with 
various observational techniques. The understanding of physi- 
cal processes in operation should also be pursued from theoret- 
ical point of view to explain the rather common occurrence of 
the warm molecular envelope, or exactly speaking, the compo- 
nent which cannot be explained in the framework of classical, 
hydrostatic photospheric models, in a wide range of late-type 
stars from Mira variables to M (super)giants and K giants. 

6. Concluding remarks 

Our simple model of the warm H2O envelope can simultane- 
ously reproduce the spectra and the angular diameters of a Ori 
and a Her obtained at 1 1 /im. For a Ori, we have also shown 
that the spectra obtained at 12 /mi as well as at 6 - 7 /im can 
be fairly reproduced by this model. The continuum-like spec- 
tra of a Ori and a Her observed in the 1 1 fim region can be 
interpreted as a result of the filling-in due to emission from the 
outer part of the warm H2O envelope. Although the observed 
featureless 1 1 //m spectra do not show any hint of the presence 
of the warm H2O envelope, it manifests itself as an increase of 
the angular diameter. Our model can reproduce the increase of 
the angular diameters observed for a Ori and a Her from the 
K band to the 1 1 /im region. For a Ori, the gas temperature, 
the H2O column density, and the radius of this H2O envelope 
were derived to be 2050 K, 2 x 10 20 cirT 2 , and 1.45 /?*, respec- 
tively. For a Her, we derived a water vapor gas temperature of 
2000 K, a column density of 7 x 10 20 cm~ 2 , and a radius of 
1.4**. 
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